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The conductance of a quantum wire containing a single magnetic barrier is studied numerically 
by means of the recursive Greens function technique. For sufficiently strong and localized barriers, 
Fano - type reflection resonances are observed close to the pinch-off regime. They are attributed to 
a magnetoelectric vortex-type quasibound state inside the magnetic barrier that interferes with an 
extended mode outside. We furthermore show that disorder can substantially modify the residual 
conductance around the pinch-off regime. 

PACS numbers: 73.23.-b,75.70.Cn 



I. INTRODUCTION 



Localized magnetic fields that are oriented perpen- 
dicular to a quantum fil m [Q , || or a quantum wire 
@, S H S 3 1, ©, EE 02, which are furthermore 
strongly localized in transport (x-) direction and homoge- 
neous in the transverse (y-) direction are known as mag- 
netic barriers (MBs). They can be realized experimen- 
tally by ferromagnetic films on top of a two-dimensional 
0, EE O [H, [H, El El or quasi one-dimensional 
electron gas residing in a semiconductor heterostructure: 
magnetizing the ferromagnetic film in x- direction results 
in a magnetic fringe field with a z-component localized 
at the edge of the film that extends along the ^-direction. 
Transport experiments on MBs in two-dimensional elec- 
tron gases show a pronounced positive magnetoresistance 
as a function of the MB amplitude El EI El El , 
which can be interpreted quantitatively in a classical pic- 
ture [TH, El , where the MB acts as a filter with a trans- 
mission probability that depends upon the anglejrf inci- 
dence of the electrons. Moreover, Peeters et al. [1] stud- 
ied the energy spectrum and the transmission properties 
of MBs with a rectangular profile in a two-dimensional 
electron gas by analytical means. This calculation pre- 
dicts resonant structures in the low energy region due to 
the presence of a virtual level Q. 

MBs in quantum wires have been the subject of several 
theoretical studies recently, which is driven by their po- 
tential ability of parametric spin filtering dUU, EHEIIj 
provided the effective g - factor is sufficiently large. MBs 
arc thus not only of fundamental interest, but also have a 
distinct potential for application in spintronics. In these 
numerical studies, resonant features in the conductance 
are frequently found, e.g in Fig . 3 of Ref. [io| . Fig. 2 of 
Ref. [IlJ or Fig. 5 of Ref. The character of these 

resonances as well as their origin has not been studied in 
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detail. This, however, is not only of fundamental inter- 
est, but also a prerequisite for possible applications. For 
example, the predicted spin polarizations can reach par- 
ticularly high values in the proximity of such resonances 

n 

Here, we use the recursive Greens function (RGF) tech- 
nique to investigate the structure and the conductance 
of a single MB that forms in a quantum wire (QWR) 
below the edge of a ferromagnetic film. The barrier 
shapes are adapted from typical experimental conditions 
[H G3 El EE M El El ■ We find that for smooth bar- 
riers with a large spatial extension, the number of trans- 
mitted modes drops stepwise and without resonances as 
the barrier amplitude increases or, correspondingly, the 
Fermi energy is reduced. As the lowest mode gets re- 
flected, the conductance of weak barriers approaches zero 
as a function of decreasing energy, a situation that we 
denote as magnetic pinch- off. For sufficiently sharp bar- 
riers, however, pronounced dips in the conductance are 
found close to the pinch-off regime. From studies of the 
local density of states (LDOS) in combination with the 
spatially resolved occupation probability densities and 
current density distributions, we conclude that the trans- 
mission zeroes originate from the interference of an ex- 
tended state with a quasi-bound vortex state that is lo- 
calized inside the MB. This leads to resonant reflection 
[20L I2lll2"2l| ar id can be regarded as a type of Fano reso- 
nance [23|, I2J]. Furthermore, we discuss the influence of 
disorder on the reflection resonances. 



II. MODEL AND CALCULATION METHOD 

Let us consider a hard- wall QWR of length L = 4 /im in 
x-direction and width w = 500 nm in y-direction, defined 
in a semiconductor heterostructure with a ferromagnetic 
film placed on its surface, Fig. [T] We use the parameters 
of GaAs, i.e. an effective electron mass of m* = 0.067m e 
in our model, where m e denotes the mass of the free 
electron. 




FIG. 1: Schematic representation of the considered two- 
terminal device, consisting of a MB created by a ferromagnetic 
film deposited on the top of a quantum wire. Also sketched is 
the square grid (period a = 5nm) used in the computation. 
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FIG. 2: The calculated conductance as a function of the Fermi 
energy Ef for MBs of different localizations (distances of the 
QWR from the sample surface, calculated for a magnetization 
of fioM = 1.2 T), and as a function of the barrier amplitude 
for Ef = 25meV (upper inset). The lower inset shows the 
shapes of the MBs present at the distances considered, for 
the film magnetization assumed in the main figure. Here, 
Ei — 22.5 /xeV denotes the ground state energy. 



The ferromagnetic film has an in-plane magnetization 
HqM in i-direction, which generates a symmetric MB 
B z (x). In- plane magnetic fields are neglected. The inho- 
mogeneous MB can be expressed as [15j 



B z (x) 



47T 



hi 



X 2 + (d+ hf 



(1) 



with h being the thickness of ferromagnetic film and d 
the distance of the QWR from the semiconductor sur- 
face, respectively. As d increases (which can be re- 
alized experimentally by preparing samples with two- 
dimensional electron gases at different locations in the 
growth (z)-direction), the localization of B z (x) is re- 
duced. A magnetization of fMoM — 1.2 T is assumed, 
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FIG. 3: The local density of states (LDOS) integrated over the 
^-direction, as a function of energy and x. (a) For d — 350 nm, 
the modes of the QWR experience a smooth diamagnetic shift 
in the barrier region, due to the formation of magnetoelectric 
subbands with spatially varying energy, (b) The LDOS for 
a strongly localized MB (d = 35 nm) shows that the smooth 
evolution of the suband energies gets interrupted in the barrier 
region, and localized states form around the center of the 
barrier. Also shown is the position of the localized state in 
relation to the conductance of the structure. 



which can be achieved experimentally by using Co [l5| 
or Dy [l8j as ferromagnetic material. It can be tuned 
by applying an external magnetic field in ^-direction 
0, UJ, 13 EE OS E3- The Fermi energy can be ad- 
justed by, e.g., a homogenous gate electrode in between 
the semiconductor surface and the ferromagnetic film. 
We consider distances of d = 15 nm, 35 nm, and 350 nm, 
which result in barriers of amplitudes B z (x = 0) = 
0.41 T, 0.28 T, 0.03 T, and full widths at half maximum 
(FWHM) of 94 nm, 148 nm, and 812 nm, respectively, see 
the lower inset in Fig. [21 Note that the integrated mag- 
netic field A = J B z (x)dx — 6.86 x 10~ 8 Tm is indepen- 
dent of d. 
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FIG. 4: Left column: gray scale plots of the sum of the prob- 
ability densities | 1 1 2 + 1 2 of the two occupied modes at 
the Fermi energy in the barrier region, shown for the energies 
marked in Fig. [2jb). The maximum of the MB is denoted 
by the dashed vertical lines, and its FWHM is 148 nm. Right 
column: the corresponding current density distributions. 



The electronic wave functions in the quantum wire ex- 
posed to the MB structure B z (x) are described by the 
effective-mass Hamiltonian 



H = Hn 



(2) 



where V c (y) is the confining potential in transverse di- 
rection which is assumed to be a hard-wall potential and 
H is the kinetic energy term. Using the Landau gauge, 
the MB can be included by choosing magnetic vector po- 
tential as A = (— B z (x)y, 0, 0). The kinetic energy can 
therefore be written as 



Ha 



2m* 



d_ 

Ox 



ieB z (x)y 



dy 2 



(3) 



In order to perform numerical computations, the com- 
putational area is discretized into a grid lattice with lat- 
tice constant a = 5 nm as shown in Fig. [TJ such that the 
continuous values x, y are denoted by discrete variables 
ma, na, respectively. Then the calculation area is con- 
nected to two ideal semi-infinite leads. The tight-binding 
Hamiltonian of the system reads 



H 



EE 



t{ C m.n c m,n+l + e 



iqw | 



tCm+l.n + H.C.} 



(4) 



where eo is the site energy which has included the ef- fects of bottom of the band and confining potential, the 
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hopping element t — h 2 /(2m*a 2 ); c\ nn and c m . n denote 
the creation and annihilation operators at the site (m, n). 
The phase factor with q = % J^ i+1 B z {x')dx' is obtained 
by using the Peierls substitution. 

In the presence of disorder, the site energy changes 
within a width A, namely 

eo — ► (o + Se Q (5) 

where the values of Seo are distributed uniformly between 
—A/2 and A/2 and A is related to the elastic mean free 
path A by [23] 

A = (6A 3 F ArVA)V2 (6) 

with Ep being the Fermi energy and the Fermi wave- 
length. 

The two-terminal conductance Git is calculated within 
the framework of the Landauer-Biittiker formalism and 
can be expressed as following [2|| 

2e 2 A 

G 2t = y ^ W 2 (7) 

0,a=l 

I 



where N is the number of propagating states in the leads, 
t/3 a is the transmission amplitude from incoming state a 
in the left lead to outgoing state /3. The transmission 
amplitudes tp a can be expressed via the total Green's 
function G of the system as tp a — ih^jv a vpG M+1 ' a , 
where G M+1 '° denotes the matrix (M + 1|G|0) with 
and M + 1 corresponding to the positions of the left and 
right leads. We calculate G M+1 '°on the basis of the re- 
cursive Green's function technique in the hybrid energy 
space formulation [U [28| . We calculate separately the 
surface Green's functions related to the left and right 
leads and the Green's function of the scattering region 
with the MB, and then link them together at the bound- 
aries. The wave function ipi for i-th slice of the region 
under consideration is calculated recursively via 

- ^ = G*°V 1( Vo + G^V^+Vi+i (8) 

with V the hopping matrix and G l0 and G" the short- 
hand notations of the Green's functions (i|G|0) and 
(i|G|i). For visualization purposes, the current density 
jmn is associated with hopping along bonds and can be 
expressed as 



jmn(r) = -i— -r{m^ n [e l9 "Vm+l,n - e i9 "Vm-l,n] + S^J^ro ,n+l ~ V'm.n-l] ~ CC.} (9) 

Zan 



where ip m n is the wave function at site (m, n), and m, n 
are the unit vectors in the longitudinal (x-) and trans- 
verse (y-) direction. 

The local density of states at the site r = (m, n) is 
related to the total Green's function in real space repre- 
sentation by the following equation [26| . 

p(r;E)=--Im[G(r,r;E)] (10) 

where Im means the imaginary part of the complex ma- 
trix. 

With this formalism, we study below a QWR with a 
hard wall confinement potential and a maximum of 4 
occupied modes outside the MB. We set the effective g- 
factor to zero, therefore the modes always have a spin 
degeneracy of 2. Via the distance d, the barrier itself is 
varied in both amplitude and localization (characterized 
by its FWHM) at fixed Fermi energy Ep. Alternatively, 
we vary Ep and keep the MB constant. In all our calcu- 
lations, the electron temperature is set to zero. 

III. RESULTS AND INTERPRETATION 

In Fig. [21 the numerical results for ballistic QWRs 
are summarized. In the upper inset, the conductance of 



a QWR with 2 occupied modes is shown as a function 
of the barrier amplitude /xqM, for distances d of 15 nm 
and 35 nm, respectively. Quantized conductance steps 
are observed, which originate from reflections of the wire 
modes at the barrier and can be regarded in close anal- 
ogy to the conductance quantization in quantum point 
contacts [29|, |30j • Above the magnetic pinch-off, the con- 
ductance shows a pronounced peak with a maximum con- 
ductance up to « 0.1e 2 /h. We note that the conductance 
between the plateau at 2e 2 /h and the peak equals zero 
within numerical accuracy. In the main figure, we show 
G as a function of Ep for the barriers given in the inset. 
For broad barriers, no structures close to the pinch-off 
regime are observed. As the localization of the MB or 
the number of occupied modes is increased, the number 
of conductance zeroes increases as well. Here, however, 
we focus on the simplest scenario where just one reso- 
nance is present. This is the reason why we study QWRs 
with at most four occupied modes. 

To shed some light on the origin of the structures in 
the conductance, we study the local density of states 
(LDOS), integrated along the y-direction, as a function 
of x (Fig. [3]). For smooth MBs (Fig. 13(a)), the MB 
generates an ^-dependent diamagnetic shift of the QWR 
modes, which are connected throughout the barrier struc- 
ture. As the Fermi energy is lowered, barriers are formed 
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for modes with subsequently lower energy, and again a 
stepwise decrease of G without resonances results. 
As the barrier is localized further, Fig. OJb), the modes 
segregate and localized states form at the center of the 
barrier, which can be regarded as remnants of the mag- 
netoelectric subbands. Qualitatively, we can understand 
this as follows. In sufficiently strong magnetic field gradi- 
ents, the magnetic phase changes abruptly in i-direction, 
and pronounced reflections occur which localize the mode 
inside the barrier. These localized states may align in 
energy with the QWR modes of a higher index and a 
resonant reflection scenario results. The formation of lo- 
calized states at the center of the MB is similar to the 
case of the double-barrier resonant tunneling (DBRT) 
structures, where transmission resonances are related to 
the existence of the quasibound states between two bar- 
riers. However, the present mechanism has a different 
phenomenology than DBRT: (i) in DBRT, the resonance 
transmission probability equals 1 for a symmetric struc- 
ture. In our system, the transmission peak is much 
smaller than 1. (ii) DBRT structures do not possess 
transmission zeroes. In our sstem, however, the dip in 
the conductance between the transmission peak and the 
first plateau equals zero within numerical accuracy. As 
the strength of the magnetic barrier increases, the height 
of the peak increases, but the transmission minimum re- 
mains at zero, (iii) For DBRT, the transmission reso- 
nances coincide with the quasi-bound states in energy 
and correspond to the poles in the complex-energy plane 
[2fj| . However, in the present case, as can be observed in 
Fig. [3] (b) (energy B), the maximum in the LDOS is at 
a different energy than the conductance peak. Rather, 
the peak position is at the low energy tail of the bound 
state. This phenomenology is the same as that one found 
in wave guides with resonators attached [2(| HH, [22j • Fol- 
lowing Shao's discussion [20| . the resonator contributes a 
phase factor A and the transmission amplitude from left 
to right can be expressed as 

Ul = t d .rl + t rs tsi/ (A- r s ) (11) 

with the transmission amplitudes t s i and t rs for being 
scattered into (from the left) and out (to the right) of the 
resonator and r s the factor from each reflection back into 
the resonator, and td, r i denotes the direct transmission 
path without a detour into the resonator. According to 
Eq. (jlip , the transmission amplitude can vanish if both 
a direct and an indirect transmission channel via the res- 
onator are present and interfere with each other. This 
behavior corresponds to a Fano resonance |23 f . which is 
also found in systems related to ours 0, l3ll. |33| . 

In the following, we argue that this scenario does ex- 
ist in a magnetic barrier under certain circumstances. 
The direct transmission channel is formed by an extended 
state which resembles an edge state in the region of high 
magnetic fields, while the bound state is a vortex state 
that forms near the center of the magnetic barrier. The 
character of such states is thus markedly different from 



snake-orbit states known from magnetic field steps with 
a change in polarity [33j. 

To further illustrate this effect, we have studied the 
spatially resolved probability density and the current 
density distribution emerging from the two occupied 
wave functions at the Fermi level close to the reflection 
resonance, see Fig. |4l where the sum of the probability 
densities \^f\ | 2 + l^l 2 of the two wave functions (belong- 
ing to the first and second energy level of the quantum 
wire) as well as the corresponding current density 
distributions are plotted as a function of the lateral 
coordinates x and y for energies A-D. In the pinch-off 
regime at energies below the transmission resonance (for 
example at energy 6.4 x E±, energy ^4), the region of high 
probability density inside the barrier is well separated 
in space from those in the leads. In addition, the 
probability of finding the electrons close to the barrier 
maximum is small. Correspondingly, the current density 
gets reflected at the flank of the barrier. As the energy 
is increased into the transmission peak (energy B), 
the region of high probability density moves along the 
x-direction into the barrier, and a significant probability 
density is found even at the barrier maximum. At the 
same time, the probability density remains asymmetric 
about the center of the QWR in y-direction. Translated 
into a current density distribution, this means that 
inside the transmission peak, a current path evolves 
where the electrons enter the barrier region at the upper 
edge of the QWR, and while a large part of the electrons 
gets rejected, a significant fraction is transmitted, via 
the lower half of the QWR, to the right hand side. 
At the reflection resonance (energy C), the region of 
high probability density is pushed even further into the 
barrier region, but at the same time develops a strongly 
symmetric shape about the center of the QWR cross 
section. This means that all the current flowing into 
the barrier gets reflected back into the QWR. At higher 
energies, the open regime is reached, like at energy D. It 
can be distinguished from the closed regime by the fact 
that a high probability density inside the barrier remains 
at only one edge of the QWR, which at the same time 
extends across the whole barrier structure in x-direction. 
This structure provides a strong transmission channel 
for the electrons, with the current flowing predominantly 
at the lower right edge across the magnetic barrier. This 
behavior is a consequence of the formation of a local 
edge state inside the magnetic barrier. 

Based on these findings, we interpret the resonances in 
the conductance as follows: close to the pinch-off regime, 
i.e. around 8E\ for d — 35 nm in our model calculations, 
there is still a small but non-vanishing direct transmis- 
sion probability through the magnetic barrier. This can 
be inferred from comparing the width of the transition 
region between the conductance plateaus 1 and 2, which 
is roughly 2.5 E\ (see Fig. [5]). The corresponding ex- 
tended state has the character on an edge state in the 
magnetic barrier. According to Eq. (jlip . it interferes 
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FIG. 5: Conductance across the magnetic barrier around the 
reflection resonance for various degrees of disorder A/Ef- 
The magnetization of the ferromagnetic film is /j,oM = 1.2 T, 
and the distance d of the QWR from the surface is 35 nm. 
Note that the parameters without disorder are identical to 
those in Fig. [2] 



have their origin in an interference between quasi-bound 
states residing inside the magnetic barrier and propagat- 
ing states of the QWR which, at the magnetic barrier, 
have the character of edge states. Even without taking 
the spin explicitly into account, it becomes clear that due 
to the resonance condition, particularly large spin polar- 
izations can be expected around the resonances, the sign 
of which should be adjustable by a small change of the 
sample parameters. More studies are necessary to an- 
alyze the details of the spin effects in such resonances. 
Furthermore, we hope that our findings stimulate exper- 
imental studies with the objective to observe this type of 
reflection resonances. 
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with the indirect transmission amplitude via the vortex- 
type bound state present in the magnetic barrier and 
generates a reflection resonance. 

We proceed by discussing the effects of disorder on 
the transmission properties. In Fig. [5l the conductance 
around the reflection resonance is shown on a logarith- 
mic scale as a function of the energy for various degrees 
of disorder. The disorder potential is modelled by a dis- 
order A of the site energy. A critical disorder energy of 
A w Ep is observed. For lower disorder, the reflection 
resonance remains basically unaffected. At larger values 
for the disorder energy, the transmission at the resonance 
becomes nonzero, while the minimum in the conductance 
shifts towards larger energies. At A/Ep « 10, the reflec- 
tion resonance is no longer a characteristic property of 
the structure, while further resonances are of compara- 
ble strength. These resonances have their origin in inter- 
ferences due to multiple reflections between impurities. 
We emphasize that while the specific disorder configu- 
ration at fixed A/Ep determines details of the energy- 
dependent conductance, it does only marginally influence 
the position, the minimum conductance or the width of 
the reflection resonance. 



IV. SUMMARY AND CONCLUSION 



The conductance of quantum wires containing a mag- 
netic barrier has been studied by the recursive Greens 
function technique. It is found that for sufficiently large 
ratios of noAl/Ep, the barrier "closes" and the trans- 
mission drop to zero. At energies close to the mag- 
netic pinch-off, reflection resonances are observed for 
sufficiently localized magnetic barriers. The resonances 
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